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Abstract. The asymptotic behavior of solutions of the three-dimensional Navier-Stokes 
equations is considered on bounded smooth domains with no-slip boundary conditions or 
on periodic domains. Asymptotic regularity conditions are presented to ensure that the 
convergence of a Leray-Hopf weak solution to its weak u-limit set - weak in the sense of 
the weak topology of the space H of square-integrable divergence-free velocity fields - are 
achieved also in the strong topology of H. In particular, if a weak w-limit set is bounded in 
the space V of velocity fields with square-integrable vorticity then the attraction to the set 
holds also in the strong topology of H. Corresponding results for the strong convergence 
towards the weak global attractor of Foias and Temam are also presented. 



1. Introduction 

The notions of limit sets and attractors (whether local or global) permeate the theory of 
dynamical systems both in finite and in infinite dimensions. In the case of infinite dimensions, 
the existence of such sets, in particular that of the global attractor, is a major issue (the 
global attractor is the minimal set for the inclusion relation which uniformly attracts all 
bounded sets of initial conditions - the global attractor contains all locally attracting sets 
and a;- limit sets). Existence results have been obtained for a number of nonlinear partial 
differential equations modeling various phenomena. 

In this note we address the celebrated system associated with the Navier-Stokes equations 
for an incompressible fluid filling a region in a three-dimensional space (3D NSE for short). 
Due in particular to the lack of a result on the global well-posedness for the 3D NSE the 
notion of attractor in this case is not settled, and the study of the asymptotic behavior of 
this system is a major challenge. 

In [7] Foias and Temam introduced a notion of weak global attractor (see the definition 
in (|3.1jl ). which is loosely speaking a global attractor for the weak topology of the natural 
phase space of square-integrable divergence-free vector fields (see also the related notion of 
trajectory attractor IH E3 ) . A notion of weak limit set (limit set for the weak topology, 
see definition (|3.3Jl )) can be similarly considered. The study of attractors and limit sets for 
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the strong topology is more delicate due to the lack of global regularity and uniqueness of 
the solutions. 

Our aim in this note is to consider weak limit sets and the weak global attractor of Foias 
and Temam and present a condition for convergence in the strong topology of the phase 
space. It is an asymptotic regularity condition. More precisely, we prove that if the weak u>- 
limit set (respectively weak global attractor) is made of points through which pass solutions 
satisfying the energy equation (with equality, not just inequality, see (|2.5|) ). then it is a strong 
tu-limit set (resp. global attractor for the strong topology). A corollary of this results with 
a simpler condition for strong convergence is that the weak u;-limit set (resp. weak global 
attractor) be bounded in the space of divergence-free velocity fields with finite enstrophy 
(i.e. square- integrable vorticity). This includes weakly attracting fixed points, weak limit 
cycles, weakly attracting quasi-periodic orbits, etc, which turn out to be strongly attracting 
provided they have bounded enstrophy. It also includes hyperbolic objects, in which case 
the weakly attracting invariant manifolds turn out to be also strongly attracting provided 
they have bounded enstrophy. 

The proof of this result is based on an idea devised by Ball (see [H EH El El ED]), exploting 
energy-type equations to prove asymptotic compactness of the trajectories. Let us consider 
limit sets for simplicity. We start with a trajectory u = u(t). This trajectory is said to be 
asymptotically compact (in a given space) if given any time sequence t n — » oo, there exists 
a convergent subsequence for {u(t n )} ra . This asymptotic compactness implies the existence 
of the corresponding tu-limit set. A similar result holds for global attractors |T2~1 1§1 [TH| FTOj . 

The idea of the energy-equation method to obtain the asymptotic compactness can be 
broken down into two steps: i) weak compactness of the sequence {u(t n )} n , and ii) norm 
convergence |u(£ n .)| — > |v | of a weakly convergente subsequence u(t nj ) — ^ v , as j — > oo. 
In uniformly convex spaces (such as Hilbert spaces), weak plus norm convergences implies 
strong convergence, hence asymptotic compactness in the strong topology. In pratice, the 
first step follows from classical a priori estimates obtained from energy-type inequalities, 
while the second one, as developed in [BE], follows from energy-type equations (the equality 
is important!). 

At this point one may be skeptical about the condition of an energy equation for the 3D 
NSE since the (Leray-Hopf) weak solutions are known to satisfy only an energy inequality. 
The crucial point that we present here is based on the simple observation that in fact the 
energy equation need only be satisfied by the limit solution associated with v . Hence, 
asymptotic regularity is all what is needed for the strong convergence of a Leray-Hopf weak 
solution towards its limit set. The precise results are given in Lemma f4. 11 Proposition 14.11 
and Theorems 14. II to 14.41 below. 

It should be clear that this idea may be adapted to yield similar results for other differential 
equations in which uniqueness and lack of regularity are troublesome, such as wave equations 
with critical nonlinearities. 
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2. Preliminaries 

We recall now some classical results which can be found for instance in fTJ EH 03 HTj . 
We consider the three-dimensional Navier-Stokes equations with either periodic or no-slip 
boundary conditions. In the periodic case, we consider the whole space M 3 , and the flow is 
assumed periodic with period Lj in each direction Xj. We define fl = I1 3 =1 (0, Li) and assume 
that the average flow on Q vanishes, i.e. 

u(x) dx = 0. 

Here, u = («!, w 2 , W3) denotes the velocity vector, and x = (xi,x 2 ,x 3 ), the space-variable. 

In the no-slip case, the flow is considered in a bounded domain Q of M 3 , with smooth 
boundary dfl, and the no-slip boundary condition on dQ is assumed, i.e. u = on dQ. 

In either periodic or no-slip case, we obtain a functional equation formulation for the 
time-dependent velocity field u = u(£) in a suitable space H: 

y + ^u + B(u,u) = f. (2.1) 

We consider the test spaces 

V = {u G C°°(Q); V ■ u = 0, u(x) is periodic with period Li in each direction Xi} , 
in the periodic case, and 

V={ueC c °°(fi) 3 ; V-u = 0}, 

in the no-slip case, and define H as the completion of V under the L 2 (Q) 3 norm (where C£°(Q) 
denotes the space of infinitely-differentiable real-valued functions with compact support on 
Q). We also consider the space V defined as the completion of V under the i? 1 (f2) 3 norm. 
We identify H with its dual and consider the dual space V, so that V C H C V. We denote 
by H w the space H endowed with its weak topology. 

We consider the inner products in H and V given respectively by 

(u ' v) = L u(x) ' v(x) dx ' ((u ' v)) = L ^ £ ' £ dx ' 

and the associated norms |u| = (u, u) 1 / 2 , ||u|| = ((u, u)) 1//2 . 

We denote by Plh the (Leray-Helmhotz) orthogonal projector in L 2 (Q) 3 onto the subspace 
H. In (12. lj) . A is the Stokes operator Au = —P^Au; B(u,v) = Plh((u ' ^) v ) i s a bilinear 
term corresponding to the inertial term; f represents the mass density of volume forces 
applied to the fluid, and we assume that f E H; and v > is the kinematic viscosity. The 
Stokes operator is a positive self-adjoint operator on H, and we denote by Ai > its first 
eigenvalue. 

A Leray-Hopf weak solution on an open time interval / = (to,ti), —oo < to < £i < oo, 
is defined as a function u = u(i) on (to,ti) with values in H and satisfying the following 
properties: 

(i) ueL TC (t ,ti;#)nL 2 oc (to,*i;V); 
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(ii) du/dteUgfahiV); 

(iii) u e C(I;Hy,) (i.e. weakly continuous); 

(iv) u satisfies the functional equation (|2.1|) almost everywhere on I = (to>*i)j 

(v) u satisfies the following energy inequality in the distribution sense on / = (i ,ti): 

^|u(*)| 2 + I /H*)|| 2 <(f,u(t)). (2.2) 

A Leray-Hopf weak solution on an interval of the form [toj^i) is defined as a Leray-Hopf 
weak solution on (to^i) which is continuous at t = to, i.e. 

(vi) u(t) — > u(t ) in H, as t — > £q . 

From now on, for notational simplicity, a weak solution will always mean a Leray-Hopf 
weak solution. For a weak solution on an arbitrary interval J, it follows that 

|u(t)| 2 < \u{t')\ 2 e-^-^ + -L|f | 2 (l - e-M^'A , ( 2 . 3 ) 

for all t in / and almost all t' in J with i' < t. The allowed times t' are the Lebesgue points of 
the function t i— > |u(i)|. In the case of a weak solution on an interval of the form [to, ii), the 
point to is a point of continuity of 1 1— > |u(t)|, hence a Lebesgue point, so that the estimate 
above is also valid for the initial time t' = t . 

Another classical estimate obtained from the energy inequality is 



|u(t)| 2 +v f ||u(s)|| 2 ds < \u(t')\ 2 + -L|f \ 2 (t - f) 
Jt' 



(2.4) 



for all t in I and almost all t' in / with t' < t, with the set of allowed times t' consisting 
again as the Lebesgue points of the function t i— > |u(£)|. 

A strong solution on an arbitrary interval I is defined as a weak solution on I satisfying 

(vii) u G C(I;V). 

Any strong solution satisfies the energy equation 

^|u(t)| 2 + i/||u(t)|| 2 = (f,u(t)) (2.5) 

in the distribution sense on its interval of definition. 

It is well established that given any initial time to and any initial condition u in H, there 
exists at least one global weak solution on [t , oo) satisfying u(t ) = u . It is also known 
that if uo belongs to V, then there exists a local strong solution, defined on some interval 
[to,t\), t\ > to, with u(io) = uo. Finally, any strong solution is unique on its interval of 
definition. The uniqueness in this case is with respect to the larger class of all weak solution, 
i.e. any weak solution with u(io) = uo must agree with the strong solution on the interval 
of definition of the latter. 
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3. Weak limit sets and the weak global attractor 

The weak global attractor, as defined in [Jj, is the set 

( there exists at least one global weak solution v = v(t) 1 

A w = < v G H; defined for all t G R which is uniformly bounded in H, > . (3.1) 
I i.e. sup igR |v(t)| < oo, and such that v(0) = v J 

Due to the energy estimate ()2.3|) and the uniform bound on the global solutions in the 
definition of A w it follows that A w is a bounded set in H: 

|v | < 4-|f|, Vv G Aw. 

It is proved in [7j that A w is weakly compact in H and that it attracts all weak solutions 
in the following sense: If u = u(t) is a weak solution on [to, oo) for some to G R, then for 
any neighborhood O of A w in the weak topology of H, there exists a time T > to such that 
u(t) G O for all t > T. Since is separable the weak topology is metrizable on bounded 
sets and the convergence above can be rewritten is terms of this metric. Finally, a certain 
invariance property holds, namely if vo belongs to A w and v = v(t), t G R, is a uniformly 
bounded global weak solution through v then v(t) G *A W for all t G R. 

Besides the pointwise attraction (attraction of individual weak solutions) of the weak 
global attractor, one can show that the attraction is in fact uniform with respect to uniformly 
bounded sets of initial condition (see |B]). More precisely, given t G R and R > 0, then for 
every neighborhood O of A w in the weak topology of H, there exists a time T > £ Q such 
that u(t) G O for all t > T and for every weak solution u on [to, oo) with sup t>to |u(t)| < R. 
Since A w is bounded in H and the weak topology of H is metrizable on bounded subsets this 
uniform attraction in the weak topology can be rewritten in terms of an associated metric. 

These properties define A w and justify its definition as the weak global attractor. They 
can also be used to characterize A w in a more classical way: 

{there exist global weak solutions u n = u n (t), n G N, defined 1 
v G H; for all t > 0, with sup tgR |u n (t)| < |f|/i/Ai, and a time > . (3.2) 
sequence {t n } n , t n > 0, t n — > oo, such that u n (t n ) — ^ v J 

Now, given an arbitrary weak solution u = u(t) on an interval of the form (to, oo) or 
[to, oo), for some to G R, we define its weak to- limit set by 

u w (u) = {v G H; 3{t„}„, t n > t , t n oo, u(t n ) -»> v } , (3.3) 

where — 51 denotes the weak convergence in H. This set is always nonempty since {u(t)} t>to 
is bounded in H (thanks to ()2.3j) ). hence weakly precompact. Since the weak topology is 

W 

metrizable on bounded subsets of H, the classical characterization u w (u) = rii>oU t > s {u(t)} 
holds, where ~ w denotes the closure in the weak topology. Hence, u w (u) is weakly compact. 
By classical dynamical system arguments one can also show that tu w (u) attracts u in the 
sense that for any weakly open set O containing cj w (u), there exists a time T > to such that 
u(t) G O for aU t > T. 



6 



RICARDO M. S. ROSA 



As for the invariance property, it is possible to show that for every vo in u w (u), there 
exists a global weak solution v = v(i), t G R, with v(0) = v and v(i) G cu w (u) for all 
i 6 1. This is achieved by passing to the limit the solutions u(t n + •) over time intervals 
[— T, T], for arbitrarily large times T. In fact, classical a priori estimates (derived from 
(Q and (Q) yield that {u(t n + •)}„> bounded on L°°(-T,T; H)(lL 2 (-T,T] V) and that 
{<9u/<9t(t n + -)} n is bounded on L 4//3 (— T, T; V), which imply precompactness in L 2 (—T; T, H) 
and C([— T, T], if w ). A diagonalization process guarantees convergence on any bounded set 
to a global weak solution defined on all R. However, due to the possible lack of uniqueness 
one cannot guarantee the invariance for every solution passing through v . 

With the above invariance property in mind, given v in uo w (x\) we defined (? u ( v o) as the 
set of all global weak solutions v = v(t), t G R, with v(0) = vo, and such that there exists 
a sequence {t n } n , t n > t , t n — > oo, with the property that u(t„ + •) converges to v in 
C([—T,T],H W ), for all T > 0. Note that this implies that Q u (^o) C u; w (u), for all v in 
w w (u). 

Finally, by classical a priori estimates (derived from (|2.3|) and (|2.4j) ) and Aubin's compact- 
ness theorem the convergence of u(t n + ■) to v holds weakly-star in L°°(— T, T; H), weakly 
in L 2 (—T,T; V), and strongly in L 2 (—T,T;H). Then, the convergence required in the defi- 
nition of Q u {v ) guarantees the uniqueness of the limit and the convergence without passing 
to further subsequences. 

A similar argument for the weak global attractor yields that for every vo in A w and every 
sequences {u n } n and {t n } n as in the characterization ()3.2j) . with u n (t n ) — ^ v , there exists 
subsequences {u n . }j and {t n . }j such that u n . {t n . + •) converges weakly-star in L°° (— T, T; H) , 
weakly in L 2 (-T, T; V), strongly in L 2 (-T, f ; H) and strongly in C(-T, T; H w ), for all T > 0, 
to a global weak solution v = v(t), with v(t) G A w , for all t G R, and with v(0) = v . 



4. Asymptotic regularity conditions for the strong convergence towards 
weak limit sets and the weak global attractor 

As mentioned in the introduction the required regularity condition is that the limit solu- 
tions satisfy the energy equation exactly. More precisely, we have the following result. 

Lemma 4.1. Let u be a weak solution defined on some interval of the form (t Q ,oo) or 
[to, oo), for some t G R. Let v G uj w (u) and let {t n } n be such that t n > t , t n — > oo, and 
u (^n) - vo weakly in H . If there exists a global weak solution v = v(t) ; i 6 1, such that 
u(t n + •) converges to v in C([—T,T],H W ), for all T > 0, and which satisfies the energy 
equation ()2.5|) on R ; then u(t n ) converges strongly in H to vq. 

Proof. In what follows, given T > we restrict ourselves to n such that t n — T > to- 

The weak solution satisfies the energy inequality ()2.2|) . By adding and subtracting the 
term i/Ai|u(t)| 2 /2, we arrive at the inequality 

^ |u(t)|2 + t 1 ^ 2 + u[u{t)]2 - (f> u(t)) ' (41) 
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where 

[u] 2 =||uf-^|u| 2 , (4.2) 

with [•] being a norm in V equivalent to || • ||. By "multiplying" (|4.1j) by appropriate 
nonnegative test functions approximating e uXt and "integrating" (e.g. use (p(t) = 9 m (t)e uXlt 
in the distribution formulation of the equation, where 9 m (t) is the continuous piecewise 
linear function joining linearly the values 9 m (—T) = 0, 9 m (—T + 1/m) = 1, 9 m (— 1/m) = 1, 
9 m {0) — 0, and let m — > oo) one finds at the limit that 

|u(t n )| 2 < |u(t n - T)|V^ T - 2 j\^ s {v[u(t n + s)] 2 - (f, u(t» + s))} ds, 

for almost all T such that t n — T > t (more precisely, it holds for all Lebesgue points T of 
the function T — > \u(t n — T)\). 

Passing to the limit as n goes to infinity and using that u(t n + •) converges to v weakly 
in L 2 (-T,T; V) we find 

limsup|u(t„)| 2 < limsup|u(t„-T)| 2 e- i/AlT -2 / e uXlS (z/[v(s)] 2 - (f, v(s))} ds. 

n^oo n— K30 J —T 

Note that we have used that [■] is an equivalent norm for V so that 

e uXlS [v(s)} 2 ds < liminf / e uXlS [u(t n + s)] 2 ds. 

-T n ^°° J-T 

As for the limit solution v it satisfies the energy equation, so that 

r° 

|v(0)| 2 = |v(-T)| 2 e-^ lT - 2 J e» XlS (z/[v(s)] 2 - (f, v(s))) ds. 

By subtracting this equality from the inequality for the limsup and using that both v(— T) 
and u(t n — T) are bounded in H independently of T and n we find 

limsup |u(i„,)| 2 - |v(0)| 2 < 2Re~ uXlT , 

for some bound R > 0. The previous passage should make clear why we need the equality 
in the energy equation only for the limit solution. Now, since T is arbitrary (except for a 
set of measure zero), we may let T — » oo to find that 

limsup |u(t n )| 2 - |v(0)| 2 < 0. 

Since v(0) = v , we obtain 

limsup |u(t n ) | < |v |. 

n^oo 

On the other hand, since u(t n ) converges weakly in H to vo we have |vo| < liminf n ->oo |u(t n )|. 
Thus lim^^oo |u(t n )| = |vo|, which together with the weak convergence implies the strong 
convergence u(t n ) — > vq in H. This concludes the proof. □ 
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Proposition 4.1. Let u be a weak solution defined on some interval of the form (to, oo) or 
[t ,oo) ; for some t G R. If for every v in u) w (u) all the global weak solutions in (? u ( v o) 
satisfy the energy equation (|2.5|) then uj w (u) attracts u in the strong topology of H. 

Proof. If this were not true we would find a time sequence {t n } n , with t n > t , t n — > oo, and 
such that {u(t n )} n does not have any subsequence converging strongly in H. Now, since 
{u(t n )} n is bounded (thanks to (|2.3Jl ) it has a weakly convergent subsequence {u(£ n .)} n ., 
with a weak limit vo belonging to u w (u). Since any solution v in (? u (vo) satisfies the condition 
of Lemma HH] it follows that {u(t n .)} n . converges strongly to v , which is a contradiction. 
Thus, co> w (u) attracts u in the strong topology of H. □ 

We now define the set 

^re g = { v o £ V'i 3<5 > and 3 a strong solution v G C((— 5, 5), V) with v(0) = v } . 

Theorem 4.1. Let u be a weak solution defined on some interval of the form (t , oo) or 
[to, oo), for some t G R. If u w (u) C V rcg! then w w (u) attracts u in the strong topology of H . 

Proof. In view of Proposition 14. II it suffices to show that any weak solution in <? u ( v o) satisfies 
the energy equation for all v in u; w (u). Note first that Q-aiyo) C V ieg since it is included in 
l) w (u), which is assumed to be included in V ieg . So any weak solution v in {? u ( v o) is included 
in V Teg , i.e. v(i) G V reg for all t G R. By the uniqueness of strong solutions among the larger 
class of weak solutions it follows that v is a global strong solution and belongs to C(R, V). 
Therefore, v satisfies the energy equation, and this completes the proof. □ 

Theorem 4.2. Let u be a weak solution defined on some interval of the form (£q,oo) or 
[<o,oo) ; for some t G R. If u> w (u) is a bounded subset of V , then u w (u) attracts u in the 
strong topology of H. 

Proof. Following the idea in the proof of Theorem 14. 11 any global weak solution v in £/ u (vo), 
v G u> w (u), is uniformly bounded in V, hence it is a global strong solution belonging to 
C(R, V), which is sufficient to apply Proposition 14. II and conclude the proof. □ 

We now present a result concerning the weak global attractor. This result is not a direct 
consequence of the previous ones simply because of the uniform attraction required (reflecting 
the fact that the weak global attractor may be larger than the union of weak w-limit sets). 

Theorem 4.3. If A w C V TCg or A w is a bounded subset of V , then A w attracts every weak 
solution in the strong topology of H , and this attraction is uniform with respect to uniformly 
bounded sets of weak solutions. More precisely, given t G R and R > 0, then for every 

e > 0, there exists a time T > t such that c?«s^(u(t), A w ) == sup V(je ^ w |u(£) — v | < e, for 
every weak solution u on [t , oo) with sup t>to |u(i)| < R. 

Proof. Suppose the result is not true. Then there exists to G R, R > 0, e > 0, a sequence 
u n of weak solutions on [t , oo) with sup t>to |u n (i)| < R, and a time sequence {t n } n , t n > t Q , 
t n —> oo, such that 

|u n (i n ) — v | > e, for all n and all v G A w . (4.3) 
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Consider the sequence v n (t) = u n (t n + t), defined for t > to — t n . By the assumption 
of uniform estimate on u n and by classical a priori estimates for the 3D NSE (derived 
from (jOD and (Q) the sequence {v n } n is bounded in L°°(-T,T; H) n L 2 (-T,T; V), and 
{dv n /dt} n is bounded in L 4 / 3 (— T, T; V), which imply precompactness in L?[— T; T, if) and 
C([— T, T], H w ). By passing to the limit as in the classical theory of existence of weak solutions 
of the 3D NSE and by using a diagonalization process we find a subsequence converging on 
any bounded interval [— T, T], T > 0, to a global weak solution v = v(i), t > 0. In particular, 
Ury(£ry) converges weakly to some element vo = v(0) in ii. 

At the limit, we retain a uniform bound for v, sup t>t() |v(£)| < R, so that v(t) belongs to 
A w for all t G R, and in particular v G A w . Since ^4 W either belongs to V Teg or is a bounded 
subset of V it follows (see the proofs of Theorems 14.11 and 14 .2|) that v is a global strong 
solution, v G C(R, V), hence v satisfies the energy equation on R. 

Now we proceed as in the proof of Lemma 14.11 The only difference being that instead of 
working with u(i„+-) we work with u n (t n + -) (we still have the inequalities holding for almost 
all T, as the countable union of zero- measure sets is still of zero measure). Then, by repeating 
the energy-equation argument with u n (t n + ■) we find that limsup^^ |u n (i n )| < |v |, so 
that in fact u n .(t n .) converges strongly to v . But this contradicts (|4.3|) since v belongs to 
A w . This completes the proof. □ 

Our last result concerns just part of the weak global attractor. It does not reduce to a 
statement about weak a»-limit sets included in the weak global attractor since it allows for 
asymptotic limits of sequences of weak solutions instead of a single solution. It applies, for 
instance, to regular (in the sense of being a strong solution) connecting orbits which are not 
necessarily cu-limit sets. 

Theorem 4.4. Suppose v = v(t), t G R ; is a global strong solution uniformly bounded in 
H (in other words, a global weak solution included in A w D V rcg ) and set v = v(0). Let 
u n = u n (t), t > 0, and {t n } n be as in the characterization (J3.2j) . with u n (t n ) — ^ v . Then, 
u n (t n ) — > Vo strongly in H. 

Proof. By assumption, v = v(t) is a strong solution, hence it satisfies the energy equation 
(J2.5j) . Then, as in the last paragraph of the proof of Theorem 14.31 we apply the energy- 
equation method to the sequence of weak solutions {u n (t n + -)} n to show that |u n (£ n )| — > |v |. 
Then, we conclude that u(t n ) vo strongly in H. (At first one may need to pass to further 
subsequences, but since the weak limit u n (t n ) — > v exists and hence is unique, the strong 
converge in H must hold for the whole sequence.) □ 
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